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We examine coherent active-active channel neutrino flavor evolution in environments where 
neutrino-neutrino forward scattering can engender large-scale collective flavor transformation. We 
introduce the concept of neutrino flavor isospin which treats neutrinos and antineutrinos on an equal 
footing, and which facilitates the analysis of neutrino systems in terms of the spin precession analogy. 
We point out a key quantity, the "total effective energy", which is conserved in several important 
regimes. Using this concept, we analyze collective neutrino and antineutrino flavor oscillation in 
the "synchronized" mode and what we term the "bi-polar" mode. We thereby are able to explain 
why large collective flavor mixing can develop on short timescales even when vacuum mixing angles 
are small in, e.g., a dense gas of initially pure v e and z/ e with an inverted neutrino mass hierarchy 
(an example of bi-polar oscillation). In the context of the spin precession analogy, we find that the 
co-rotating frame provides insights into more general systems, where either the synchronized or bi- 
polar mode could arise. For example, we use the co-rotating frame to demonstrate how large flavor 
mixing in the bi-polar mode can occur in the presence of a large and dominant matter background. 
We use the adiabatic condition to derive a simple criterion for determining whether the synchronized 
or bi-polar mode will occur. Based on this criterion we predict that neutrinos and antineutrinos 
emitted from a proto-neutron star in a core-collapse supernova event can experience synchronized 
and bi-polar flavor transformations in sequence before conventional Mikhyev-Smirnov-Wolfenstein 
flavor evolution takes over. This certainly will affect the analyses of future supernova neutrino sig- 
nals, and might affect the treatment of shock re-heating rates and nucleosynthesis depending on the 
depth at which collective transformation arises. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

In both the early universe and in core-collapse super- 
novae, neutrinos and antineutrinos can dominate ener- 
getics and can be instrumental in setting compositions 
(i.e., the neutron-to-proton ratio). However, the way 
these particles couple to matter in these environments 
frequently is flavor specific. Whenever there are differ- 
ences in the number fluxes or energy distribution func- 
tions among the active neutrino species (y e , v e , v^, v^, v r 
and v T ), flavor mixing and conversion can be important 
[1-8]. 

In turn, the flavor conversion process becomes compli- 
cated and nonlinear in environments with large effective 
neutrino and/or antineutrino number densities [1, 3, 5, 9- 
11]. In these circumstances neutrino- neutrino forward 
scattering can become an important determinant of the 
way in which neutrinos and antineutrinos oscillate among 
flavor states. 

Two of the three vacuum mixing angles for the active 
neutrinos are now measured. The third angle (#13) is 
constrained by experiments and is limited to values such 
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that sin 2 26*13 ;$ 0.1 (see, e.g., Ref. [12] for a review). In 
addition, the differences of the squares of the neutrino 
mass eigenvalues are now measured, though the abso- 
lute masses and, therefore, the neutrino mass hierarchy 
remains unknown. 

Both the solar and atmospheric neutrino mass- 
squared differences are small, so small in fact that con- 
ventional matter-driven Mikhyev-Smirnov-Wolfenstein 
(MSW) evolution [13-15] would suggest that neutrino 
and/or antineutrino flavor conversion occurs only far out 
in the supernova envelope. On the other hand, it has 
been shown that plausible conditions of neutrino flux 
in both the early shock re-heating epoch and the later 
neutrino-driven wind, r-process epoch, could provide the 
necessary condition for neutrino-neutrino forward scat- 
tering induced large-scale flavor conversion deep in the 
supernova environment [16]. 

The treatment of the flavor evolution of supernova neu- 
trinos remains a complicated problem, and the exact so- 
lution to this problem may only be revealed by full self- 
consistent numerical simulations. However, physical in- 
sights still can be gained by studying somewhat simpli- 
fied models of the realistic environments. For example, 
one source of complication is that there are three active 
flavors of neutrinos in play. As the measured vacuum 
mass-squared difference for atmospheric neutrino oscilla- 
tions (i5m 2 tm ~ 3 x 1CP 3 eV 2 ) is much larger than that for 
solar neutrino oscillations (5m% ~ 8x 10 -5 eV 2 ), the gen- 
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oral problem of three-neutrino mixing in many cases may 
be reduced to two separate cases of two-neutrino mixing, 
each involving v e (y e ) and some linear combination of 
and v T [Dfjt and v T ). This reduction allows the possibil- 
ity of visualizing the neutrino flavor transformation as 
the rotation of a "polarization vector" in a three dimen- 
sional flavor space [17]. Different notations have been 
developed around this concept (see, e.g., Refs. [3, 18]). 
However, none of these notations fully exhibits the sym- 
metry of particles and anti-particles in the SU(2) group 
that governs the 2x2 flavor transformation. 

The equations of motion (e.o.m.) of a neutrino "po- 
larization vector" is similar to those of a magnetic spin 
precessing around magnetic fields. One naturally expects 
that some collective behaviors may exist in dense neu- 
trino gases just as for magnetic spins in crystals. Indeed, 
it was observed in numerical simulations that neutrinos 
with different energies in a dense gas act as if they have 
the same vacuum oscillation frequency [10]. This collec- 
tive behavior was later explained by drawing analogy to 
atomic spin-orbit coupling in external fields and termed 
"synchronization" [6] . 

Another, more puzzling, type of collective flavor trans- 
formation, the "bi-polar" mode, has been observed in 
the numerical simulations of a dense gas of initially pure 
v e and D e [19]. This type of collective flavor transfor- 
mation usually occurs on timescales much shorter than 
those of vacuum oscillations. Although the analytical so- 
lutions to some simple examples of "bi-polar" systems 
have been found [20, 21], many aspects of these bi-polar 
systems still remain to be understood. In particular, it 
seems counter-intuitive that, even for a small mixing an- 
gle, large flavor mixing occurs in both the neutrino and 
antineutrino sectors in a dense gas initially consisting of 
pure v e and v e for an inverted mass hierarchy. 

Both synchronized and bi-polar flavor transformation 
were discovered in the numerical simulations aimed at 
the early universe environment. It has been shown that 
synchronized oscillation can also occur in the supernova 
environment [22]. However, it is not clear if supernova 
neutrinos can also have bi-polar flavor transformation. 
If supernova neutrinos can have collective synchronized 
and/or bi-polar oscillations, the questions are then where 
these collective oscillations would occur and how neutrino 
energy spectra would be modified. 

In this paper we try to answer the above questions. 
In Sec. II we will give the general equations govern- 
ing the mixing of two neutrino flavors in the frequently 
used forms and introduce the notation of neutrino flavor 
isospin, which treats neutrinos and antineutrinos on an 
equal footing. We will also point out a key quantity, the 
"total effective energy", in analogy to the total energy 
of magnetic spin systems, which is conserved in some in- 
teresting cases. In Sec. Ill and Sec. IV we will analyze 
the synchronized and bi-polar neutrino systems using the 
same framework in each case. We will first describe and 
explain the main features of these collective modes using 
the concept of total effective energy. We then generalize 



these analyses by employing "co-rotating frames" . We 
will derive the criteria for the occurrence of these collec- 
tive modes, and discuss the effects of an ordinary matter 
background. In Sec. V we will outline the regions in 
supernovae where the neutrino mixing is dominated by 
the synchronized, bi-polar and conventional MSW fla- 
vor transformations. We will also describe the typical 
neutrino mixing scenarios expected with different matter 
density profiles. In Sec. VI we will summarize our new 
findings and give our conclusions. 

II. GENERAL EQUATIONS GOVERNING 
NEUTRINO FLAVOR TRANSFORMATION 

We consider the mixing of two neutrino flavor eigen- 
states, say \v e ) and \v T ), which are linear combinations 
of the vacuum mass eigenstates \u\) and |^) with eigen- 
values mi and m^, respectively: 

\v e ) = cos# v |z/i) + sm9 v \v2), (la) 
\v T ) = - sm6 v \ui) + cos# v |^ 2 ), (lb) 

where 9 V is the vacuum mixing angle. We take 9 V < 7r/4 
and refer to dm 2 = m 2 — m\ > as the normal mass 
hierarchy and dm 2 < as the inverted mass hierarchy. 
When a neutrino with energy E v propagates in matter, 
the evolution of its wavefunction in the flavor basis 

^(:::) < 2 > 

is governed by a Schrodinger-like equation 

ij t ^ = (H v +H e )^, (3) 

where a Vc and a Ur are the amplitudes for the neutrino to 
be in \v e ) and \u T ) at time t, respectively. (This equa- 
tion is "Schrodinger-like" because, unlike the Schrodingcr 
equation, we arc here concerned with flavor evolution at 
fixed energy and with relativistic leptons.) The vacuum 
mass contribution H v to the propagation Hamiltonian in 
the flavor basis is 

<5m2 / _ cos 29 v sin 26> v \ . . 

~ AE V \ sin20 v cos26» v ) ' W 

and the contribution H e due to forward scattering on 
electrons in the same basis is 

«.=!(; .!)• < 5 » 

where A = ^/2G-pn e with n e being the net electron num- 
ber density. Eq. (3) also applies to the antineutrino wave- 
function 




3 



if A in H e is replaced by —A. 

When a large number of neutrinos and antineutrinos 
propagate through the same region of matter, their for- 
ward scattering on each other makes another contribu- 
tion to the propagation Hamiltonian for each particle. 
For the ith neutrino, this contribution is [1, 3, 5, 9, 23] 



J 



D 

b: 



D, 



(7) 



where 



B = V2Gf ^(1 - COS©y){n^[(p^)ee - {Pv,])tt] - np.j[{pp.j)ee - (Pp^tt]}, (8a) 

j 

B eT = 2\/2G F y^(l - cos 8jj )[n„,j(p„,j) eT - n 9t j(p V j) eT \. (8b) 
j 

I 



In the above equations, 0^ is the angle between the prop- 
agation directions of the ith neutrino and the jth neu- 
trino or antineutrino, and n v j (n„j) and p v j (pp,j) are 
the number density and single-particle flavor-basis den- 
sity matrix of the jth neutrino (antineutrino), respec- 
tively. Specifically, 



and 



Pv 



a* a VT 



Kl 5 



a v ,a„ 



(9a) 



(9b) 



where we have adopted the convention for the density ma- 
trix of an antineutrino in Ref . [3] . The neutrino- neutrino 
forward scattering contribution for an antineutrino can 
be obtained by making the substitution B — > —B and 



B e 



-B* in H 



The single-particle density matrices in Eqs. (9a) and 
(9b) can be written in the form 



\(l + P-<r), 



where P is the polarization vector in the three- 
dimensional (Euclidean) flavor space and a represents 
the Pauli matrices. Explicitly, the polarization vectors 
in column form are 



P„ = 



2Re« e a„ T ) 
2Im(a* e a lv ) 
\a v \ 2 - \a v J 2 



(11) 



and 



2Re(a Pc ap T ) 
2Im(ap e ap ) 
\a,\ 2 -\aj\ 2 



(12) 



By straightforward algebra, it can be shown that the 
Schrodinger-like equation 



(13) 



(10) and a similar equation for an antineutrino lead to [3] 




V2G F n e | 



j 

- \FlG-pne | | -V2Gf ^(1 - cos ®ij){n v ,jP v ,j - n p jP p j) 



, (14a) 
(.14b) 



The three real components of the polarization vector which is irrelevant for flavor transformation. Therefore, 
contain the same information as the two complex am- Eqs. (14a) and (14b) are equivalent to the Schrodingcr- 
plitudes of the wavefunction except for an overall phase like equations. Eqs. (14a) and (14b) appear to suggest 
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a geometric picture of precessing polarization vectors. 
This picture has been discussed quite extensively in the 
literature (see, e.g., [6, 24, 25]) and shown to be es- 
pecially helpful in understanding flavor transformation 
when neutrino self-interaction (i.e., neutrino-neutrino 
forward scattering) is important. To facilitate the use 
of this picture, we briefly discuss the physics behind it 
and introduce some notations. 

For simplicity, we first consider only the contributions 
H v and H e to the propagation Hamiltonian H for a neu- 
trino. In this case, we can write 



= ■ (a*vHv 



H e ), 



where 



H v 

He 

if 



Sm 2 

— sin 29 v + e z cos 28 v , 



IV2G 



(15) 



(16) 

(17) 
(18) 



with e4 and e* being the unit vectors in the x- and z- 
directions in the flavor basis, respectively. Eq. (15) takes 
the form of the interaction between the "magnetic mo- 
ment" fi — 7s of a spin-i particle and an external "mag- 
netic field" H = H cff / 7 with 



H 



off 



He 



(19) 



Here 7 is the "gyromagnctic ratio" and can be chosen 
arbitrarily. Classically, the spin s would experience a 
torque t = /xxH = sx H cff and its e.o.m. would be 
given by the angular momentum theorem: 



d 

dt S 



s x H 



off 



(20) 



By Ehrenfest's theorem, the quantum mechanical de- 
scription of a system has the same form as the classical 
e.o.m. provided that all physical observables are replaced 
by the expectation values of their quantum mechanical 
operators. In the present case, if we replace s in Eq. (20) 

by 



(21) 



then neutrino flavor transformation governed by H 
H v + H e can be described quantum mechanically by 



dt 



= s„ x H 



eff 



(22) 



expectation value s„ of this operator as the NFIS and 
use it instead of the polarization vector P„ to describe 
neutrino flavor transformation. The z-componcnt of a 
NFIS s„ is of special importance as it determines the 
probability for the corresponding neutrino to be in \v e )'- 



3 vz — S ^ 



■ e, = 



= a,. - 



(23) 



Therefore, for a neutrino, s i vz = 1/2, —1/2 and corre- 
spond to \u e ), \v T ) and a maximally mixed state, respec- 
tively. 

Adiabatic MSW flavor conversion has a simple expla- 
nation in this "magnetic spin" analogy. For illustrative 
purposes we assume Sm 2 > and ff, C 1. As a y e prop- 
agates from a region with large matter density, e.g., the 
core of the sun, to a region of very little ordinary mat- 
ter, H cff changes its direction from ~ —e z to Hy ~ e*. 
If the density of electrons n e changes only slowly along 
the way (adiabatic process), H eS also changes slowly, and 
the NFIS s„ corresponding to the neutrino is always anti- 
aligned with H off . Therefore, the neutrino originally in 
the v e eigenstatc (s„ = e z /2) is now mostly in the u T 
eigenstate (s„ ~ —e z /2). 

It is useful to illustrate the criterion for adiabadicity 
of this process in the "magnetic spin" analogy. First, 
we note that the probabilities for a neutrino to be in 
instantaneous mass eigenstates vl (light) and vyt (heavy) 
are 



1 + cos 26» 

2 ' 
cos26» 



(24a) 
(24b) 



respectively, where 28 is the angle between the directions 
of s„ and e™ = H ofi /|H eff |, and e™ are the unit vectors 
for the instantaneous mass basis. In the MSW picture, 
is the instantaneous matter mixing angle. In an adiabatic 
process, \a VL \ 2 and |a„ H | 2 are constant, and so is 9. Using 
Eq. (22) we have 

~(co S 20) = ^.e 7) = Sv .^. (25) 

On a timescale St > 2n/\H cfi \, s v has rotated by at 
least one cycle around H off . If H cff changes its direc- 
tion only by a small angle 8<p = |de™/dt|5t <C 2tt during 
St, then s„ in Eq. (25) averages to (s„ • e™)e™. Noting 
that ef ■ (def/dt) = (l/2)d(|e™ | 2 )/dt = 0, one can sec 
that the angle 6 is unchanged in this process. There- 
fore, the criterion for a MSW flavor transformation to be 
adiabatic is 



which is the same as Eq. (14a) in the absence of neutrino 
self-interaction. Clearly, the operator a/2 in Eqs. (15) 
and (21) represents a fictitious spin in the neutrino flavor 
space, which may be appropriately called the neutrino 
flavor isospin (NFIS). The flavor eigenstates \v e ) and \v T ) 
correspond to the up and down eigenstates, respectively, 
of the z-component of er/2. We will loosely refer to the 



dt 2 



H cff x H 
|H cff l 2 



cir 



< |H 



off 



(26) 



which is equivalent to saying that the rate of change of 
the direction of the "magnetic field" H cff is much smaller 
than the rotating rate of the "magnetic spin" s„ around 
H eff . 
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The full version of Eq. (14a) can be obtained by ex- 
tending the Hamiltonian in Eq. (4) to include 



(27) 



where 

My - = -2V2G F (1 - cosGy). (28) 
We define the NFIS for an antincutrino as 1 



S;, = -- 



Pp 

2 ' 



(29) 



so that the terms related to neutrinos and antineutrinos 
appear symmetrically in H. vv .i- The probability for an 
antineutrino to be in \v e ) is determined from 



f _ -f 



an - an 



1 , 

2 ~ |ac 



(30) 



For an antineutrino, s^ z = 1/2, —1/2 and correspond 
to \P T ), \v e ) and a maximally mixed state, respectively. 

Now Eqs. (14a) and (14b) can be rewritten in terms of 
the NFIS's in a more compact way 



— s 4 = s t x ( /iv,*H v + H e + ^2 I':,"--,*, ) > ( 31 ) 



with the understanding that 



_ J Sm 2 / (2E Ut i) for a neutrino, 



-dm /(2Ep.i) for an antineutrino, 



(32) 



and that the sum runs over both neutrinos and antineu- 
trinos. We also define a total effective energy (density) £ 
for a system of neutrinos and antineutrinos that interact 
with a matter background as well as among themselves 
through forward scattering: 

£ = ~ ^ ] TT-i/.i^i ' — (JiijTli/jTlvjSi • Sj, (33) 

* hj 

where 



H® ff = ^iv,iHv + H e 



(34) 



We note that this effective energy should not be confused 
with the physical energies of neutrinos and antineutrinos. 
It can be shown from Eq. (31) that £ is constant if n e and 



1 The two fundamental representations 2 and 2 of the SU(2) group 
generated by the Pauli matrices are equivalent. These represen- 
tations are related to each other by the transformation a y , and 
tpi? = (Tytpi? transforms in exactly the same way as does ip v under 
rotation. Defining sa = i/>t(<x/2)i/>n , one naturally obtains the 
minus sign in Eq. (29). 



all the n v i's and also constant. The concept of 

the total effective energy will prove useful in understand- 
ing collective flavor transformation in a dense neutrino 
gas. 

In the early universe the neutrino gas is isotropic, and 



"ij —> [i-v = -2v / 2Gf- 



(35) 



For illustrative purposes we will assume this isotropy con- 
dition in most of what follows. We will discuss the ef- 
fects of the anisotropic supernova neutrino distributions 
in Sec. V. 



III. SYNCHRONIZED FLAVOR 
TRANSFORMATION 

In a dense neutrino gas NFIS's are coupled to each 
other through self-interaction and may exhibit collec- 
tive behaviors. As discovered in the numerical simula- 
tions of Ref. [10], neutrinos with different energies in a 
dense gas act as if they are oscillating with the same 
frequency. This collective behavior was referred to as 
"synchronized" flavor oscillations in the literature and 
explained in Ref. [6] by drawing analogy to atomic spin- 
orbit coupling in external magnetic fields. In this section 
we will first review the characteristics of a simple syn- 
chronized NFIS system from the perspective of the con- 
servation of the total energy £ of the NFIS system. We 
will then extend the discussion to more general synchro- 
nized NFIS systems using the concept of a "co-rotating 
frame" and demonstrate the criteria for a NFIS system 
to be in the synchronized mode. We will show that the 
stability of a synchronized system is secured by the con- 
servation of the total effective energy. In the last part of 
the section we will look into the problem of synchronized 
flavor transformation in the presence of ordinary matter, 
which is relevant for the supernova environment. 



A. A Simple Example of Synchronized Flavor 
Transfer mat ion 

We start with a simple case of a uniform and isotropic 
neutrino gas with no matter background (n e = 0). The 
gas initially consists of pure neutrinos with a finite energy 
range corresponding to |/Uv,i| < |Mv,»|max, and all the 
n v /s stay constant. The e.o.m. of a single NFIS is 



dt 



Si — Si x (^v,iH v + Mi,S), 



where 



(36) 



(37) 



is the total NFIS (density) of the gas. (The NFIS den- 
sity for an individual "spin" is just n u jSj.) Summing 
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Eq. (36) over all neutrinos, we obtain 



B. General Synchronized Systems 



-^-S = ^ HV,i n v,i S i X H V- 



(38) 



Following the discussion at the end of the preceding sec- 
tion, the evolution of the individual (sj) and the total (S) 
NFIS obeys conservation of the total effective energy 

£ = — /iV,i n iAi s i • Hy 

i 

— ^ ^ f-iijTiv.iTifs jSi • Sj (39a) 

i,3 

= - ^2l^v,in v ,iSi ■ H v - ^S 2 - ( 39b ) 
An interesting limit is 

,% | max- 

(40) 

Noting that each Sj has a magnitude of 1/2 and Hy has 
a magnitude of unity, we see that 



£~-^S 2 ~ const. 



(41) 



in the above limit. Therefore, a gas with a large initial 
total NFIS S evolves in such a way that it roughly main- 
tains the magnitude of its S. For such a gas, Eq. (36) 
reduces to 



dt 



&i — f-t>i/&i X S, 



(42) 



which means that each precesses around the total NFIS 
with a fixed common (angular) frequency 



ui v = \fJ, v S\. 



(43) 



Eq. (38) shows that S evolves on a timescale > 
27r/|/Uv,i|max- Consequently, over a period St satisfying 



2tt r 2?r 

< St < : , 

|MV,i | max 



(44) 



Si averages out to be (s^ • S)S/S 2 and Eq. (38) effectively 
becomes 



— S ~ w sync S x H v , 



where 



^sync — (Mv) = ^ 



(45) 



(46) 



It can be shown from Eqs. (38), (41), and (42) that 
const. Therefore, S processes around Hy with 



^sync 



a fixed frequency ui sync while the individual Sj's precess 
around S with a fixed common frequency u u . This collec- 
tive behavior of a dense neutrino gas is usually referred 
to as synchronized flavor oscillations [6]. 



Synchronization can occur not only in dense neutrino 
gases but also in dense antincutrino gases and gases in- 
cluding both neutrinos and antineutrinos. Noting that 
the NFIS's for neutrinos and antineutrinos essentially 
only differ by the signs in /J,y/s [see Eq. (32)], one can 
repeat the same arguments in Sec. Ill A for these more 
generalized cases. Instead of doing so, we want to pro- 
ceed from a new perspective, which demonstrates some 
of the benefits of the NFIS notation. 

We consider a reference frame rotating with an angular 
velocity of — flHy. In this co-rotating frame, Eqs. (36) 
and (38) take the form 

Si = Si x (/* v ,iH v + n„S), (47a) 
S = /jy,in„,iSi x H v , (47b) 



where s, (S) and §j (S) are Sj (S) and its time derivative 
in terms of their x-, y-, and z-components in the co- 
rotating frame, and 



(48) 



It is clear that one can set fly of a NFIS to any value by 
choosing an appropriate co-rotating frame, and a NFIS 
for an antineutrino in the lab frame becomes a neutrino in 
some co-rotating frame. For example, the NFIS in the lab 
frame with s = — e l z /2 and ^y = —5m 2 /2E corresponds 
to a D e with energy E. In a co-rotating frame with = 
-Sm 2 /E the NFIS has s = -e f z /2 and fly = Sm 2 /2E, 
which corresponds to a v T with energy E. Therefore, the 
NFIS notation really treats neutrinos and antineutrinos 
on an equal footing. 

Because S and S are the same vector in two different 
frames, the synchronization of the NFIS's in one frame 
means the synchronization in any frame. Consequently, 
synchronization can occur in dense antineutrino gases 
and gases of both neutrinos and antineutrinos just as 
it can occur in pure neutrino gases as long as Eq. (40) is 
satisfied in some co-rotating frame. 

As we have seen, l/Lty^lmax is not uniquely determined 
and can have different values in different co-rotating 
frames. However, we note that the relative spread of the 
individual values of the fiy/s of the NFIS's is an intrin- 
sic property of a NFIS system and is co-rotating frame 
invariant. For a co-rotating frame with 



one has 



where 



(MV,i)min + (MV,i) r 



|Av,i| max — A^iy, 



a _ (MV,i)max ~ (MV,i) r 

A Mv = ■ o 



(49) 



(50) 



(51) 
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measures the spread of the /ly.i's in the NFIS system. 
Synchronization can be obtained if 



|M„S| > A/Zy. 



(52) 



When applying this condition to astrophysical envi- 
ronments such as the early universe and supernovae, 
we must consider the meaning of A^y as neutrinos 
in these environments formally have an infinite energy 
range. One interesting scenario is where the distribution 
of NFIS density as a function of fiy has a single dominant 
peak. An example is the neutronization burst in a core- 
collapse supernova event where the neutrinos emitted are 
dominantly v e with a Fcrmi-Dirac-like energy distribu- 
tion f Ve (E). For this case a natural estimate of A/xy 
is the half- width of the distribution function f Us {^v), 
where /i/ c (mv) is obtained from f Vc {E) using the relation 
Hv(E) = 8m 2 /2E. 

Another interesting scenario is where the distribution 
of NFIS density function of /iy has two domi- 

nant peaks. An example of this scenario is the Kelvin- 
Helmholtz cooling phase of a proto-neutron star in a core- 
collapse supernova event where the neutrinos emitted are 
mostly (in number) v e and v e . For this scenario one can 
take 



A^iy 



(5m 2 <5m 2 
+ 



2E V 2Ep c 



(53) 



where E Ve and i^ Pe are the peak energies of the v e and 
v e energy spectra, respectively. 

For more complicated scenarios, the criterion to ob- 
tain synchronized flavor oscillations can be compared to 
the criterion for an adiabatic MSW flavor conversion. If 
a NFIS system has been tested to be in a synchronized 
mode using the analyses in Sec. Ill A in some co-rotating 
frame, each individual NFIS s should precess around the 
total NFIS S with a fixed angle. This is the same "track- 
ing" behavior as in the adiabatic MSW flavor transfor- 
mation process discussed in Sec. II except that S now 
takes the place of H cff in Eq. (26). Because S slowly ro- 
tates around Hy with angular frequency Lu sync = (/ty), 
the adiabatic criterion yields 



|/i„S|>|<£v>Hsin£|, 



(54) 



where £ is the angle between the directions of S and Hy. 
Eq. (54) provides a necessary condition for synchroniza- 
tion. Practically one may use 



|^s|>|( M v)| 



(55) 



as the criterion for synchronization, where (/iy) is eval- 
uated using Eq. (46) with all the relevant neutrino and 
antineutrino energy distributions. 

We now make some comments on the stability of the 
synchronized mode. Because neutrinos with different en- 
ergies have different vacuum oscillation frequencies, one 
may think that the NFIS's will develop relative phases 



and that the resulting destructive interference will break 
the synchronization, i.e., reducing S to approximately 0. 
Indeed, using Eq. (38) one can see that 



dt 



S 2 = 2S • S 



(56a) 



^ n v ^n v j{nY,i - Mv,j)(sj x s t ) ■ H v (56b) 



is generally not zero, and therefore, |S| varies with time. 
However, Eq. (36), from which Eq. (38) is derived, can 
be used to show that the total effective energy £ is 
conserved and the total NFIS S roughly maintains con- 
stant magnitude if the n v /s do not vary with time, [see 
Eq. (41)]. In this case, destructive interference stemming 
from the relative phases of different NFIS's cannot com- 
pletely destroy synchronized flavor oscillations. On the 
other hand, if |S| ~ initially, no significant synchro- 
nization of NFIS's can occur spontaneously. This result 
is in accord with the lengthy study in Ref. [26]. 



C. Synchronized Flavor Transformation with a 
Matter Background 

We now discuss the effects of a matter background on 
synchronized flavor transformation in dense gases of neu- 
trinos and/or antineutrinos. The relevant evolution equa- 
tions are 



dt 
d_ 

dt 



Si X (/iy.jHv + H e + ^„S), 



(57a) 



S = X! Mv^jSj x H v + S x H e . (57b) 



First we assume a fixed matter background with net elec- 
tron number density n e . For high n e , corresponding to 
|H e | > |(/U V )|, Eqs. (57a) and (57b) reduce to 



_d 

dt 
d_ 

dt 



Si ~ Si x (H e + n„S), 

S ~ S X H e . 



(58a) 
(58b) 



The above equations correspond to perfectly synchro- 
nized flavor oscillations: in a frame rotating with an an- 
gular velocity of H e , the total NFIS stays fixed and 
the individual NFIS's precess around it with a common 
frequency |/i„S|. However, for neutrinos and antineutri- 
nos initially in pure flavor eigenstates, s* and S start out 
aligned or anti-aligned with H e = — e^y^Gpne- There- 
fore, the above perfect synchronized flavor oscillations 
reduce to a trivial case where all Sj's remain in the initial 
state (i.e., all neutrinos stay in their initial flavor states). 
This trivial case is of no interest to us and will not be 
discussed further. 

For |H e | ~ |(mv)I an d |/z^S| » |(/Uy)|, the discussion 
is similar to the case with no matter background. All Sj's 
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precess around S with a frequency |/x„S| and Eq. (57b) 
becomes 



-S~Sx((/x v >H v + H e ). 
at 



(59) 



Therefore, the total NFIS of the gas precesses around the 
effective field H cfi = (/xv)Hv + H e and behaves just as 
does a single NFIS with s = S/(2|S|) and /x v = (/xv) 
in the same matter background [see Eq. (22)]. For the 
cases with Sm 2 > and (/xv) > or with Sm 2 < 
and (/xv ) < 0, this representative NFIS corresponds to a 
neutrino with energy 



E = 

-^sync — 



5m 2 



2(/x v ) 



(60) 



For the other cases, this representative NFIS corresponds 
to an antineutrino with energy E sync . For an initially 
pure v e neutrino gas, E~ y \ Q is simply the neutrino energy 
distribution- averaged value of E~ x : 



E- 1 

-^sync 



USE) 
E 



dE, 



(61) 



where f Ve (E) is the energy distribution of v e . For more 
general cases, E sync is evaluated using Eqs. (46) and (60) 
with all the relevant neutrino and antineutrino energy 
distributions. 

The above discussion can be extended to the case of 
a slowly varying matter background in a straightforward 
manner. We note that this is again an adiabatic process 
as discussed in Sec. II except that S takes the place of 
s„ this time. The angle between S and H cfi is there- 
fore constant. A gas of initially dominantly v e with 
|/x„S| 3> |(/xv)| acts just like a single neutrino with en- 
ergy -Esync propagating in this matter background. For 
a normal mass hierarchy (dm 2 > 0), there may be an 
MSW resonance that can enhance flavor transformation. 
In contrast, no MSW resonance exists and flavor transfor- 
mation is suppressed by the matter effect for an inverted 
mass hierarchy (Sm 2 < 0). 

Obviously, for a neutrino and/or antineutrino gas with 
|/x^S| <C |(/xv)|, there is no synchronized flavor transfor- 
mation. 



IV. BI-POLAR FLAVOR TRANSFORMATION 

The astrophysical environments where neutrino flavor 
transformation is of interest do not always provide condi- 
tions which are favorable for synchronization. For a neu- 
trino gas to be in the synchronized mode, the neutrinos 
have to be prepared in such a way that the corresponding 
NFIS's are strongly aligned in one direction. There are 
important regimes where this does not occur. 

For example, consider the 2x2 mixing channels v e ^ 
v T and v e ^ D T in the late-time, shocked region above the 
proto- neutron star. By definition, Sj = e^/2 for a v e or v T 
and Si = —e z /2 for a D e or v T . Therefore v e , 9 e , v T , and 



v T form two NFIS blocks pointing in opposite directions 
when they leave the neutrino sphere. The subsequent 
behavior of these neutrinos is interesting. We show be- 
low that under the right conditions large-scale collective 
"swapping" of flavors v e ^ v T and D e ^ D T can occur 
in a mode in which the NFIS blocks remain more or less 
oppositely-directed. This is an example of the bi-polar 
mode. 

In Ref. [19], numerical simulations of a homogeneous, 
dense neutrino-antineutrino gas in the absence of a mat- 
ter background showed that the flavor "swapping" in the 
bi-polar mode occurred at a higher frequency than would 
vacuum oscillations. Ref. [20] gave an analytical solution 
to a simple bi-polar system, a gas initially consisting of 
equal numbers of mono-energetic v e and D e , for a normal 
mass hierarchy. Ref. [21] generalized the solution to a 
gas of unequal numbers of v e and v e with different en- 
ergies, again for a normal mass hierarchy scenario, and 
found that the system exhibits bimodal features (dual 
frequencies) . 

In this section we again adopt a physical, analytical 
approach and use a simple example to illustrate neutrino 
flavor transformation in bi-polar systems from the energy 
conservation perspective. For the first time, we explain 
why large flavor mixing can develop in some bi-polar sys- 
tems even with a small mixing angle. We will then extend 
the discussion to more general bi-polar systems using the 
co-rotating frame, and discuss how bimodal features can 
appear in such systems. We will propose criteria under 
which a NFIS system can be in the bi-polar mode, and 
show that a bi-polar system is at least semi-stable. We 
will conclude this section with some discussion on the 
effects of the matter background on bi-polar flavor trans- 
formation. 



A. A Simple Example of Bi-Polar Flavor 
Transformation 

We start with a simple bi-polar system initially con- 
sisting of mono-energetic v e and v e with an equal number 
density n v , which form two NFIS blocks Si(0) = n„s„ e = 
e l z n u /2 and S 2 (0) = n„Sp e = —e z n u /2. This system is 
uniform and isotropic and has no matter background. 
The evolution of Si and S2 is governed by [see Eq. (31)] 



dt 
d 
dt 



Si = Si x (/x v ,iH v + /^ s 2), 



(62a) 



S 2 = S 2 x (mv,2Hv + j*„Si), (62b) 



where /tv,i = — A*v,2 = Mv- With the definition of 

S+ = Si + S 2 and S_ = Si-S 2 , (63) 

we find 

d 

d? 
d_ 
dt' 



-S + = /x v S_ x H v , (64a) 
: S . = /x v S+ x H v + /^S_ x S+. (64b) 
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The initial conditions are 
0, 



S+(0) 
S_(0) 



(65a) 



n v e z = n„(e^ s'm29 v + e v z cos2# v ), (65b) 



where and are the unit vectors in the x- and z- 
directions, respectively, in the vacuum mass basis (ej = 
Hy). Using these conditions and Eq. (64), we can show 
that 



= s_.e; = o. 



(66) 



In other words, S+ can only move parallel to while S_ 
is confined to move in the plane defined by and e v z . 

The evolution of S + and S_ obeys conservation of the 
total effective energy 



s/ 2 



£ — — fly_\S\ • Hy — fly ,2^2 ' Hy 
h) 

= — ^yS_ • Hy — ^-S^_ = — ^yn„cos20 v , 
which gives 

|S_| cost? = n„ cos 20 v - -^-Si 

2/x v 



(67a) 
(67b) 

(68) 



with § being the angle between S_ and e* z . Further, it 
can be shown from Eq. (63) that 



S 2 + + S 2 _=nl 
Combining Eqs. (68) and (69), we obtain 
cos20 v fx v n v ( 1 



COS"!? = 



where 



2/A 



|S_|/n„. 



(69) 



(70) 



(71) 



Noting that fx v = —2^/2Gf < and s_ < 1, we see 
that for a normal mass hierarchy {fxy > 0), cos2(? v < 
cost? < 1 and S_ is constrained to oscillate around e v z 
with — 28 v < t? < 29 v . For S v Cl, the system stays close 
to the initial state and there is very little flavor mixing. 

The situation for an inverted mass hierarchy (fiy < 
0) is more complicated. We proceed by first rewriting 
Eq. (70) as 



cos i? = ^cos 28 v — 
= cos20 v 

where 



fi v n v 

2 fly 



1 

+ 

s_ 
1 



2mv 



2^ v 



cos 2B V 



s_ (72a) 
(,72b) 

(73) 




FIG. 1: The solution of si for a simple bi-polar system in the 
vacuum mass basis for a normal (solid line) and an inverted 
(dashed line) mass hierarchy, respectively. The solution of S2 
can be obtained from that of Si using Eq. (75). 



is a positive characteristic neutrino number density for 
fiy < 0. The evolution of the bi-polar system under 
consideration falls into the following three categories de- 
pending on the parameter n u /n"\ For n v jn c J l < 1/2, 
d cos t?/ds_ < and S_ is constrained to oscillate around 
with — 29 v < t? < 29 v just as in the case of a normal 
mass hierarchy. In this case, the difference between fiy^\ 
and fiy ^ is too large for the two NFIS blocks to main- 
tain strong correlation, and neutrinos and antineutrinos 
oscillate as two separate sectors. For 1/2 < n„/n™ < 1, 
dcost?/ds„ is positive for s_ ~ 1 but becomes and 
then negative for smaller s_, and the maximum value 
t?max corresponding to dcost?/ds_ = is given by 



cos t? max = 2 cos 26 v * 



„cri 



(74) 



In this case, t? first increases from 28 v for the initial state 
to t? max as s_ decreases and then decreases to as s_ 
further decreases to its minimum value. Subsequently the 
motion of S_ is mirrored in the other half of the plane 
defined by and e v z . Note that for each complete cycle 
S_ reaches the position at t? = 29 v (—26 v ) twice but with 
S- = 1 and a smaller value, respectively. For n^/n" 1 > 
1, Eq. (72) shows that dcost?/ds_ is always positive, 
— 1 < cost? < cos2# v and S_ oscillates around with 
29 v < t? < 2tt-26» v (note that for n„/n™ = 1, S_ shrinks 
to at t? = 7r/2 and 37r/2, therefore appearing to skip 
the range tt/2 < t? < 3tt/2). In the limit n„/n™ > 1, s_ 
stays ~ 1 as S_ rotates in the plane defined by and 
e v z . This is particularly interesting because the two NFIS 
blocks remain anti-aligned and can completely reverse 
their initial directions, which means that full conversion 
of the initial v e and v e occurs even for 9 V <C 1 . 

The above results on the evolution of S_ and S+ can be 
adapted easily to describe the evolution of the individual 
NFIS's si and S2 of an initial v e and v e , respectively. 
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FIG. 2: (Color online) The evolution of the z-component of si in the flavor basis for the simple bi-polar system with (a) 
a normal and (b) an inverted mass hierarchy. The solid lines are for |/iv,i|/|/ii/^| = 1/10, and the dashed lines are for 
I/Am \/\fJ-i>n v — 1/40. The vacuum mixing angle 8 V is taken to be 0.1. 



Using Eqs. (63) and (66), we find 



From the above equations and Eq. (75), we obtain 



sin$ 



2n„ 



'12/ 



'lz 



'22/ 



>2z 



2n v 
2n v 



S- COS1? 



(75a) 
(75b) 
(75c) 



where for example, s\ x is the x-component of Si in the 
vacuum mass basis. As a visual illustration, we show in 
Fig. 1 the evolution of Si in the vacuum mass basis for 
the case of a normal mass hierarchy (solid curve) and for 
the case of an inverted mass hierarchy with n„/n™ > 1 
(dashed curve). The trajectory of S! for each case marks 
the intersection between the parabolic surface represent- 
ing 



/A' 



cos 2ft 



J + ix v n v {sl y ) 2 = 0, 



(76) 



which is equivalent to Eq. (68), and the spherical surface 
representing 



(sl x ? + {sl v ) 2 + {sl z f 



(77) 



which follows from the fixed mag nitude of = (1/2) 2 . 
The evolution of S2 can be obtained from that of si based 
on Eq. (75). 

Of course, the flavor evolution of an initial v e is de- 
scribed most directly by the z-component s\ z of s\ in the 
flavor basis. The unit vectors in this basis are related to 
those in the vacuum mass basis as 



eZ cos 2ft, — eY sin 2ft, 



eY sin 2ft, + eY cos 2ft, 



(78a) 
(78b) 

(78c) 



sl x sin 20 v + s\ z cos 26 v 



cos(i?-20 v ). (79) 



As S_ evolves very little from the initial state for 8 V <C 1 
in the case of a normal mass hierarchy and in the case of 
an inverted mass hierarchy with n^/n"' 1 < 1/2, there is 
little flavor evolution in these cases. In contrast, for the 
case of an inverted mass hierarchy with n u /n c J l 3> 1, we 
have s_ ~ 1 and 29 v < < 27r — 2ft,, so an initial v e can 
be converted essentially fully into a v T even for S v < 1. 

Taking ft, = 0.1 and \fj, v \n u /\^ v \ = 10 (n^/n™ = 5.1 
for < 0), we show the time evolution of s\ z as the solid 
lines in Figs. 2(a) and 2(b) for a normal and an inverted 
mass hierarchy, respectively. The cases with the same 8 V 
but \ii v \n v /\liv\ = 40 (n„/n™ = 20.4 for /a/ < 0) are 
shown as the dashed lines. [In order to show the small 
evolution in the case of a normal mass hierarchy, we have 
greatly expanded the vertical scale in Fig. 2(a).] 

We note that the period of vacuum oscillations in these 
numerical examples is T vac = 27t/|/a/|- This is longer 
than the bi-polar oscillation periods Tbi shown in Fig. 2. 
In addition, the bi-polar oscillation periods decrease by 
a factor of 2 when the neutrino density n„ is increased 
by a factor of 4. These observations can be understood 
from Eq. (64), even without an outright solution of this 
equation. In the limit |/Xj/|n„/|/A/| ^ 1> the second term 
on the right hand side of Eq. (64b) dominates, and S_ 
simply rotates around S + with roughly a constant mag- 
nitude n v and frequency 



The average value of |S+| 
estimated from Eq. (64a): 



l/M 



(80) 



in the above equation can be 



l/A/S- 



l/A/K- 



(81) 
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TABLE I: The conditions for a dense neutrino gas starting 
as two groups of mono-energetic neutrino species with equal 
number to develop large flavor mixing in the small mixing 
angle scenario. Combinations of neutrino species other than 
those shown here will not develop large flavor mixing in this 
case. 







v T —v T 


V e -V T 




8m 


> 


Never 


Always 


E Ve > E Vt 


Ev T < Ev e 


Sm 2 


< 


Always 


Never 


E Uc < E Vt 


Ey T > E De 



Combining Eqs. (80) and (81) we obtain 

Tbi ~ 1 , ■ (82) 

This simple dimensional analysis agrees with the exact 
expression for the bi-polar period in the normal mass hi- 
erarchy case [20] . Therefore, for a large neutrino density 
n v >• \fiy/fi u \, the bi-polar oscillation period Tu is much 
smaller than T vac and decreases as <~ 

B. General Bi-Polar Systems 

We next look at a slightly more complicated neutrino- 
antineutrino system. In particular, we consider a system 
which is the same as that discussed in Sec. IV A except 
that v e and v e have different energies. We again define 
S + and S_ as in Eq. (63). Using Eq. (62), we find 

-^S+ = S+xH+ + S_xH_, (83a) 
-^S_ = S_ x H+ + S+ x H + /i„S_ x S+,(83b) 
where 

H ± = " v ^ v ' a H v . (84) 

When viewed in the reference frame rotating with angular 
velocity — H + , the e.o.m. of S + and S_ derived from 
Eq. (83) are exactly the same as that in Eq. (64) and 
must, therefore, have the same solution. When viewed in 
the lab frame, this NFIS system not only demonstrates 
the bi-polar oscillation as discussed in Sec. IV A, but also 
rotates around Hy at the same time. We regard this 
kind of flavor transformation as also being of bi-polar 
type. Note that this bi-polar system has two intrinsic 
periods, i.e., Tu and 27r/|H + |. This bimodal feature of 
the neutrino-antineutrino system was first discussed in 
Ref. [21]. 

Note that the above arguments employing co-rotating 
frames apply not only to systems consisting of v e -v e , but 
also to systems of v T -v Tl v e ~v T or v T -v e . Because v e -v e 
systems can develop large flavor mixing in the case of a 
small mixing angle and an inverted mass hierarchy, the 



other systems also can exhibit the same phenomenon as 
long as 

A*v,i < Mv,2, (85) 

where /xy,i and nv,2 are the vacuum coupling coefficients 
of the "spin-up" and "spin-down" NFIS's, respectively. 
For convenience, we have listed these conditions in Ta- 
ble I. 

From the simple examples discussed above we can infer 
a general description of a system possessing bi-polar os- 
cillations: a system composed of two groups of neutrinos 
and/or antineutrinos of roughly equal numbers, where 
the corresponding NFIS's point in two roughly opposite 
directions and have different characteristic values of /iy. 

In Fig. 3 (the solid lines) we illustrate two non-ideal bi- 
polar systems. These examples consist of gases of initially 
pure v e and v e with n Ve ^ n Pe and E Ve ^ Ep c . One of 
the examples [Fig. 3(a)] illustrates neutrino mixing with 
a large mixing angle and Sm 2 ~ SttIq, and the other 
[Fig. 3(b)] uses a small mixing angle and Sm 2 ~ —5m 2 tm . 

As the difference between the densities of the two NFIS 
blocks becomes larger and larger, one of the NFIS blocks 
will eventually dominate the other, and the system will 
become synchronized rather than bi-polar. This can been 
seen from Eq. (83). For simplicity, we work in the frame 
rotating with angular velocity H + where the c.o.m. of 
S± take the same form as Eq. (64). We note that a key 
characteristic of any bi-polar system is a configuration 
in which a large and near constant magnitude S_ vector 
rotates about S+ [Eq. (64b)]. In this configuration, S + 
typically has a small, variable magnitude. If one of the 
NFIS blocks dominates, it is possible that |/i„S+| will be 

bigger than |S+ x S + |/|S+| 2 , and the adiabatic condition 
can be satisfied. (Note that the rate of change of S + is 
bounded by the intrinsic frequency of the bi-polar oscil- 
lation.) If the adiabatic condition is satisfied, then S_ 
will precess rapidly around S+, with a constant relative 
angle between them. In this case, S_ will average out to 
be (S_ • S+)S + /|S + | 2 . At the same time, S + will have 
roughly constant magnitude and will rotate around Hy 
with a angular frequency [Eq. (64a)] 

* + - (^1^) (86) 

Clearly, the dominant oscillation behavior of the NFIS 
system is the slow rotation around Hy and the synchro- 
nized mode obtains in this case. Indeed, one can explic- 
itly show that Cj+ is the synchronization frequency ui sync 
of the system in the co-rotating frame. 

From the above arguments one can see that the condi- 
tion for a bi-polar system to degrade into a synchronized 
mode is the same as that for S__ to adiabatically precess 
around S+. Therefore, the criterion for a NFIS system 
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FIG. 3: (Color online) The evolution of the z-component of si (in the flavor basis) for representative bi-polar systems. These 
systems consist of mono-energetic neutrinos in initially pure and v e flavor states. The mixing parameters 8 V and Sm 2 are 
0.56 and 8 x 10 -5 eV 2 , respectively, for the calculations in panel (a), and are 0.1 and —3 x 10 -3 eV 2 , respectively, for panel (b). 
The energies of v e and v e are taken to be HMeV and 16MeV, respectively. The (effective) number densities of v a and are 
10 28 cm" 3 and 6.9 x 10 27 cm -3 , respectively. The solid, dashed, and dot-dashed lines are plotted with electron number density 
taken to be 0, 10 26 , and 10 29 cm -3 , respectively. 



to be bi-polar is that 

I/VS+I < |S | 5 X | S+I = |^s y „c| • |sin^|, (87) 

l s +l 

where £ is the angle between the directions of S + and 
Hy. We note that this criterion is exactly the opposite of 
the synchronization criterion given in Eq. (54). In many 
cases, the criterion for bi-polar flavor transformation can 
be expressed as 

|/i,S|<|( M v)|. (88) 

We now comment on the stability of the bi-polar mode. 
Realistic systems of interest usually consist of neutri- 
nos and/or antineutrinos with continuous energy dis- 
tributions. Because neutrinos (antineutrinos) of differ- 
ent energies have different vacuum oscillation frequen- 
cies, one might suspect that the bi-polar mode eventu- 
ally collapses as a result of destructive interference. We 
have argued above (Sec. IIIB) that the conservation of 
total effective energy £ essentially guarantees the sta- 
bility of the synchronized mode. This conclusion does 
not extend directly to the bi-polar mode. However, this 
energy conservation condition does provide some shield- 
ing of the two-oppositely-directed-NFIS-block configura- 
tion against rapid destructive interference-driven break- 
down. This is because the effective energies of the two 
NFIS blocks (-^Sf/2 ~ and -^Sl/2 ~ 

— /z y n^ 2 /8) and the interaction energy of these blocks 
(— /x„Si • S2 ~ Hui^i/, 1^,2/4) sum to almost zero, and be- 
cause each of these ingredient energies are large in mag- 
nitude and not easily altered. If the bi-polar configura- 
tion is ever to break down, the two NFIS blocks must 



disassemble simultaneously in a symmetric way in order 
to conserve the total effective energy. As a result, the 
bi-polar mode is at least semi-stable. It has been ob- 
served that the two-oppositely-directed-NFIS-block con- 
figuration is roughly maintained in example numerical 
simulations [18, 27]. 

C. Bi-Polar Flavor Transformation with a Matter 
Background 

In Fig. 3 we present results of numerical solutions to 
the e.o.m. for two NFIS blocks [Eq. (31)]. In this figure, 
we show examples (the dashed and dot-dashed lines) of 
flavor oscillations in gases of initially pure v e and P e in the 
presence of various matter backgrounds. It can be seen 
that flavor transformation is suppressed for the scenario 
with a large mixing angle and a normal mass hierarchy 
if n e > n v . On the other hand, large flavor mixing still 
occurs in these systems with small mixing angles and an 
inverted mass hierarchy even if the electron density dom- 
inates, although the flavor oscillation period is somewhat 
longer than in the case with no matter background. 

This phenomenon can be understood qualitatively by 
using the concept of co-rotating frames. In the presence 
of a matter background, the motions of S± still are gov- 
erned by Eq. (83) except that H + is in this case defined 
as 

H + = ^^ V - 2 H V + H, (89) 

We decompose H_ into two components: H j_ and 
H_ 11. These vectors are perpendicular and parallel to 
H + , respectively. In the reference frame rotating with 



13 




H 




FIG. 4: A bi-polar system in the presence of a large electron 
background, (a) H- can be decomposed into two compo- 
nents, H_,x and HLj, which are perpendicular and parallel 
to H+, respectively, (b) Component H x will rotate very 
rapidly and will have little net effect in the frame rotating 
with angular velocity — H+. In the same co-rotating frame, 
H_ i || is static. This component drives the bi-polar motion of 
the system. 



angular velocity — H + , we have 



H_,± = -H_ j_ x H + , 
H_ = 0. 



(90a) 
(90b) 



This configuration is illustrated in Fig. 4. If n e is very 
large, H_ j_ will rotate very rapidly and the NFIS's are 
not able to follow it. In this limit, H_ x will have on av- 
erage negligible influence on the overall evolution of the 
system, at least so long as Si and S2 are not aligned 



with H 



■II 



Note that this scenario is similar to the 



simple small mixing angle bi-polar example discussed in 
Sec. IV A. However, one difference is that here H_,|| 
takes the place of /iyHy in the simple case. Bi-polar 
systems with matter backgrounds, therefore, behave sim- 
ilarly to those without. 

Consider the simple v e -v e system discussed in 
Sec. IV A but now with a large matter background. The 
two NFIS blocks formed by neutrinos and antineutri- 
nos are initially aligned or anti-aligned with H e = H + . 
The field component H_,j_ will perturb the system to 
break this alignment. The perturbation by H_,j_ is not 
enough to cause Si and S 2 to deviate much from their 
original directions. This is because H_x has negligi- 
ble net effects once the angle between Si (S2) and H_ 11 
is significant. However, the configuration with a mis- 
alignment of Si and S 2 relative to H_,|| is just like 
the initial configuration of the simple v e -v e system dis- 
cussed in Sec. IV A. Therefore, for the normal mass 
hierarchy (5m 2 > 0), Si and S 2 will oscillate around 
*- f - II. II 



el = 



H 



-■II 



H 



and no significant fla- 
vor mixing occurs. For the inverted mass hierarchy 



-■II 



(5m 2 < 0), Si and S 2 could completely swap their di- 
rections and large flavor mixing will result. 

Dynamically, the configuration with S_ aligned with 
H_ 11 is "stable" because H_m acts as a restoring "force" 
which prevents S_ from becoming significantly mis- 
aligned with it. This is why neutrino- antineutrino gases 
of initially pure v e and v e do not have large flavor mixing 
if Sm 2 > and 5 V < 1. On the other hand, the con- 
figuration with S_ anti-aligned with H_,|| is "unstable" 
because H_,|| acts as a "force" which drives S_ toward 
the ("stable") position of alignment with H_ ^ and then 
drives it back to the original ( "unstable" ) position. This 
is why neutrino-antineutrino gases of initially pure v e and 
v e can have large flavor mixing if Sm 2 < and V <1. 



V. COLLECTIVE NEUTRINO FLAVOR 
TRANSFORMATION IN SUPERNOVAE 

To treat neutrino and antineutrino flavor transforma- 
tion in a core-collapse supernova event, we must account 
for nonuniformity and anisotropy in the neutrino den- 
sity distribution. At a radius r which is larger than the 
neutrino sphere radius i?„, the neutrino number density 
distribution is 



d 2 n v 



L v f v (E v ) 





, o<e<e , 

otherwise, 



(91) 



where L v is the neutrino (energy) luminosity, f v (E v ) is 
the normalized neutrino energy distribution, (EJ) is the 
average neutrino energy, df2„ is the differential solid angle 
around the radial direction with being the polar angle, 
and 



cosGq = VI - (Ru/rf 



(92) 



In general, flavor evolution of neutrinos traveling in 
different directions above the neutrino sphere will be dif- 
ferent due to the anisotropy of the neutrino density dis- 
tribution. For a qualitative discussion, we will assume 
the "single-angle approximation" (see, e.g., Ref. [16]) 
that the flavor evolution history of a radially propagat- 
ing neutrino is representative of all neutrinos. Under this 
approximation, 



ff„(E v )&E v 



where 



.off 



/< 



1-cosG) 



&E v dSl v 



dn u dE v 



AnRl(E v ) 



1 - y/1 - (R v /rf 



(93) 

(94a) 
(94b) 



We will comment on the validity of the single-angle ap- 
proximation at the end of the section. 
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When they leave the neutrino sphere, the neutrinos 
and antincutrinos form two oppositely directed NFIS 
blocks. For illustrative purposes we assume that v e and 
v e dominate the neutrino species emitted from the proto- 
neutron star. We will also assume that v e and v e have 
the same luminosity L v , and that the vacuum coupling 
coefficients of the two corresponding NFIS blocks are 

6m 2 dm 2 
|Uv,i — 777^; — 7 and ^ v ,2 — —- 



(95) 



respectively. 

Wc define a dimcnsionless quantity 

\5m 2 \/2(E„) 

\n„\nf 
\Sm 2 \TrRl 

~ 3.6 x 1CT 6 



1-y/l- (R„/r) 2 
| <5m 2 1 



3 x lO- 3 eV 2 
10 51 erg/s' 



10 km 



1-y/l- {R v /rf 



(96a) 
(96b) 

(96c) 



This quantity gives a measure of the inverse of the num- 
ber density of cither neutrino species. Using Eqs. (88), 
(94b), (95) and (96) we find that the rough boundary 
condition for supernova neutrinos to transition from the 
synchronized mode to the bi-polar mode is 



where the dimensionless quantity 

_ ((E Ve )-(E Pe )) 2 



2((E l , c ) 2 + (E [ ,J 2 )- 



(97) 



(98) 



measures the disparity between the energy spectra of v e 
and v e . If R v is much smaller than the boundary radius 
res (Bi-polar Starting) of the two collective modes, we 
can estimate 



res — (51 km) 



\8m 2 



-1/4 



3 x 10" 3 eV 2 

V2 , „ N 1/4 



10 51 erg/s 



10 km 



(o.oi) 



1/4 



(99) 



As neutrinos propagate away from the proto-neutron 
star, the local neutrino density decreases. Beyond some 
radius the neutrino density is so low that the collectivity 
of neutrino flavor transformation breaks down and neu- 
trinos undergo conventional MSW flavor evolution. This 
occurs if 



A//- 



vi (2) 



> 



fj,„S 



1(2)1 



(100) 



where A^ V i(2) is the half- width of the distribution 
fv e (v e ){n-v) (see the discussion in Sec. IIIB). We esti- 
mate that 



AMvi(2) ~ 



\5m 2 



2(£ M p e )) 



2 A£ "e{?«)' 



(101) 



where AE v jp e ^ is the half-width of the v e (y e ) energy 
spectra. Using Eqs. (94b), (96), (100) and (101) we can 
obtain a condition for where collectivity of neutrino flavor 
oscillations will break down: 



> (E v 



2AE V 



(102) 



If R v is much smaller than the boundary radius where 
collectivity breaks down, tbe (Bi-polar Ending), we can 
show that 



r BE (193 km) ( 



\5m 2 



V3 x 10- 3 eV 2 



-1/4 



10 51 erg/s 



R v y/AEr/WY 1 " 



10 km 7 V 0.25 



1/4 



(103) 



Taking \5m 2 \ = 3 x 10" 3 eV 2 , R v = 10km, L v = 
10 51 erg/s and AE V /(E V ) = 1/4, we can calculate the 
boundary radius between the two collective modes and 
the radius of the boundary separating the collective 
modes from the regime where conventional MSW evolu- 
tion dominates. These boundaries are shown in Fig. 5. It 
is clear that supernova neutrinos are in the synchronized 
mode near the proto-neutron star (region I in Fig. 5), but 
could experience bi-polar flavor transformation at a mod- 
erate distance (region II) . It is only in the region far from 
the proto-neutron star that neutrinos will undergo con- 
ventional MSW transformation (region III). This is very 
different from the solar neutrino oscillation case, where 
neutrinos experience only MSW flavor transformation. 
We note that there are actually no sharp boundaries be- 
tween these flavor transformation regions. There will be 
neutrinos and antineutrinos with many different energies 
in any region above the proto-neutron star. As a result, 
a particular region in general could host superpositions 
of various neutrino and antineutrino oscillation modes. 
Therefore, regions I, II and III should be understood 
as where synchronized, bi-polar, and conventional MSW 
type flavor transformations, respectively, dominate. 

In broad brush, the mixing parameters (0 V and 8m 2 ), 
the neutrino and antineutrino energy spectra and lumi- 
nosities, and r and R v are the principal determinants of 
the dominant oscillation mode at a particular location. 
Of course, the actual detailed form of flavor oscillation at 
any point is also affected by the matter density and elec- 
tron fraction. For example, if v e number flux dominates 
over that for v e and 5m 2 > 0, neutrinos and antineutrinos 
in the synchronized mode will evolve as if they were one 
neutrino with energy E sync [Eq. (60)]. As a result, neu- 
trinos and antineutrinos in this case will experience flavor 
conversion at radius ~ ^MSw(-E sy nc), where rMSw(-Ev) is 
the MSW resonance radius for a neutrino with energy 
E v . This radius is determined by the standard MSW 
resonance condition, 



Sm 2 
2EZ 



cos26> v = v / 2G , Fn e (rMSw(^)). 



(104) 
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FIG. 5: The regions where the neutrino oscillations in the 
supernova environment are dominated by synchronized (I), bi- 
polar (II) and conventional MSW (III) flavor evolution. In the 
calculations for this figure we have taken \Sm 2 | = 3 x 10~ 3 eV 2 , 
R v = 10 km, L v = 10 51 erg/s and AE V /{E„) = 0.25. 

Governed by the density run in the supernova envelope, 
neutrinos and antincutrinos may experience the following 
collective flavor mixing scenarios: 

• If the proto-neutron star has a very "thick" enve- 
lope, n e is very large throughout regions I and II 
and rMSw(-Esync) ^ ^be- In this case no significant 
flavor conversion will occur when neutrinos are in 
region I. For Sm 2 > flavor conversion is also sup- 
pressed in region II. For Sm 2 < 0, however, large 
flavor mixing can occur in region II. 

• If the envelope of the proto-neutron star is "thin" , 
n e is very small in regions II and III, and 
^MSW^sync) ;$ »"bs- F° r bm 2 > 0, almost com- 
plete flavor conversion {v e — > u T and v e — ► D T ) oc- 
curs around radius rMSw(^sync) in region I. Enter- 
ing region II, neutrinos and antincutrinos are now 
dominantly v T and v T , and large flavor mixing will 
occur (see Table I). For 8m 2 < 0, flavor conversion 
is suppressed in region I, and large flavor mixing 
can occur in region II. 

• If the proto-neutron star has an envelope of a mod- 
crate thickness, n e is large in region I and part of 
region II and r B s < r M sw(E syn c) < ^be • In this 
case, flavor mixing is always suppressed in region I. 
For Sm 2 < 0, large flavor mixing can occur in re- 
gion II. The exact flavor oscillation form is not clear 
for the 5m 2 > case in region II. However, some 
resonance- like behavior around radius rMSw(-E'sync) 
could be expected. 

The typical collective oscillations described above are 
idealized. In particular, we have assumed that the gradi- 
ent of n e is not so large that the adiabatic condition is vi- 



olated. If the adiabatic condition is violated, other inter- 
esting phenomena may occur. For example, the NFIS's 
of neutrinos and antincutrinos may be kicked into a con- 
figuration where they are aligned or anti-aligned with e { x 
at some instant. This corresponds to a maximally mixed 
state. If this occurs in region I and the adiabatic con- 
dition holds from this point on, the NFIS's will rotate 
collectively around H cfi = (/lyJHy + H e . This nearly 
maximal mixing will last until radius res or rMSw(-^sync)) 
whichever is smaller. This is the Background Dominant 
Solution described in Ref. [16]. 

Our analysis of collective flavor transformation as- 
sumed isotropy of the neutrino gases. This is appropriate 
for the early universe scenario. In the supernova environ- 
ment the neutrino gas is not isotropic. The anisotropy 
of the supernova neutrinos has two major effects. One 
effect is that the neutrinos scattering at some particular 
point have travelled different distances from the neutrino 
sphere before they interact. We note that the propaga- 
tion distances along various trajectories are most differ- 
ent near the neutrino sphere. Close to the proto-neutron 
star, the electron density is so large that n e 3> n° ff , and 
H e breaks the correlation of the NFIS's on different tra- 
jectories. As a result, the NFIS's on different trajectories 
develop relative phases. This effect, however, does not 
compromise our analysis because neutrinos and antineu- 
trinos are essentially kept in their flavor eigenstates by 
H e , and the effects of destructive interference are small. 
After neutrinos propagate away from the proto-neutron 
star, the distance difference between any two trajectories 
becomes small. 

The other effect of the anisotropy of supernova neutri- 
nos is that the neutrino-neutrino forward scattering po- 
tential, and therefore /i^ , depends on the angle between 
the directions of the neutrino momenta [Eq. (28)]. As a 
result, the effective total NFIS 

Sf e^^^^Sj (105) 

3 

is different for NFIS's on different trajectories, and one 
cannot define a universal total NFIS S in the original 
isotropic sense. However, we note that the NFIS's on 
different trajectories are still strongly coupled as a result 
of large neutrino density. This is why collective neu- 
trino flavor oscillations can arise in the first place. Al- 
though the exact neutrino oscillation behavior can only 
be shown by the numerical simulations which treat the 
trajectory issue self-consistently, we expect that qualita- 
tively similar flavor oscillations, e.g., large flavor mixing 
in the Sm 2 < and 6 V <C 1 scenario, may occur in the 
real supernova environment. 

VI. CONCLUSIONS 

We have introduced a notation for neutrino flavor 
isospin which explicitly exhibits symmetry between flavor 
transformation of neutrinos and antincutrinos. We have 
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pointed out a key quantity in dense gases of neutrinos 
and/or antineutrinos, the total effective energy, which is 
conserved in some interesting cases. Using the conserva- 
tion of the total effective energy, we have proved the sta- 
bility of synchronized flavor transformation in a simple 
and intuitive fashion. We have also demonstrated how 
co-rotating frames can be useful in analyzing collective 
oscillation in more general cases. 

With the concept of total effective energy we have for 
the first time explained why large flavor mixing occurs 
for a dense gas of initially pure v e and v e with a small 
mixing angle and an inverted mass hierarchy We have 
estimated the oscillation periods of bi-polar systems us- 
ing simple dimensional analysis. Additionally, we have 
studied more complicated and more general bi-polar sys- 
tems by using co-rotating frames. We have also for the 
first time demonstrated that a dense gas initially con- 
sisting of pure v e and v e with an inverted mass hierarchy 
can develop large flavor mixing, even in the presence of 
a dominant matter background. 

We have derived a convenient criterion for determining 
whether the synchronized or bi-polar type of collective os- 
cillations may arise in a dense neutrino and/or antineu- 
trino gas. Based on this criterion, we have estimated the 
regions where various modes of flavor oscillation may oc- 
cur in the supernova environment. We have found that 
neutrinos emitted from the proto-neutron star in a corc- 
collapse supernova event generally experience synchro- 
nized and bi-polar flavor transformations in sequence be- 
fore the conventional MSW flavor transformation takes 
over. We have also described the typical flavor oscilla- 
tion behaviors according to different density runs in the 
supernova envelope. 

Although our analysis of neutrino flavor transforma- 
tion in the supernova environment is based on crude es- 
timates, it does suggest a picture of neutrino flavor trans- 
formations dramatically different from that in the solar 
case. In particular, because of the large neutrino lumi- 



nosities, both synchronized and bi-polar types of collec- 
tive flavor transformations are involved in the supernova 
scenario. 

To go beyond this work we would need to drop a num- 
ber of the approximations made here and go over to a 
detailed numerical simulation starting from realistic con- 
ditions of neutrino and antineutrino luminosity and spec- 
tral distribution. Chief among the requirements of a de- 
tailed numerical model would be a self-consistent treat- 
ment of flavor evolution on different trajectories from the 
neutrino sphere. This could be especially important for 
regions near the neutrino sphere. Even when such de- 
tailed numerical simulations are accomplished, the col- 
lective behavior of neutrino flavors will remain a com- 
plicated phenomenon. This is where our simple physical 
pictures may be most useful: delineating the expected 
qualitative behavior of the self-interacting neutrino sys- 
tem in various supernova conditions. 

In any case, our results have probably overturned some 
of the existing paradigms related to the supernova neu- 
trino flavor oscillation problem. Among these, the analy- 
ses of future supernova neutrino signals are certainly af- 
fected, because most if not all current analyses are based 
on the assumption that the conventional MSW transfor- 
mation is valid throughout the supernova environment. 
Depending how deep the collective large-scale flavor mix- 
ing of neutrinos and antineutrinos may occur, the treat- 
ment of shock re-heating and nucleosynthesis might also 
be affected. 
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